function [] = displayEstimates(runId)

%%% Add paths
addpath('../general_funcs');

%%% Get input / output paths
config = getConfig_pairAdoptions(runId);
runFolder  = config.runFolder;

load(sprintf('%s/results.mat', runFolder), 'params_hat', 'params_ses', 'modelFits');

%%% Load data to get paramNames
load(sprintf('%s/runData.mat', runFolder));
paramNames = get_param_names_pairAdoptions(data);
NumObs = size(data.Y(:),1) / 2014;
clear data;

% Store LL and NumParams
LL = modelFits.LL;
AIC = modelFits.AIC;
BIC = modelFits.BIC;
NumParams = length(paramNames);

% Make table of estimates and display it
M1 = reportSubmodelResults(paramNames, params_hat, params_ses, modelFits, {}, true);
clear params_hat params_ses modelFits;


LL
AIC
BIC
NumParams
NumObs

M1b = cell2table(formatNumbersWith1000commas([NumParams;NumObs]));
M1b.Properties.RowNames = {'NumParams', 'NumObs'};
M1b.Properties.VariableNames = {'estimateStrings'};
M1 = [M1; M1b];
disp(M1)

writetable(M1, sprintf('%s/results2.csv', config.runFolder), 'WriteRowNames', true);

end
